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In the last few years, there has been an Increasing interest in 
the theory of differential equations with retarded arguments or systems 
in which the rate of change of a system may depend upon its past history. 

This is partially due to the fact that such equations arise in a natural 
manner in certain types of control problems. Much of the recent literature 
has been devoted to the extension of known results for ordinary differential 
equations to differential equations with retarded arguments. The present 
paper is another step in this direction. 

More specifically, we shall indicate in what manner a particular 
form of the method of averaging of Krylov-Bogoliubov-Mitropolski-Diliberto 
can be extended to differential equations with hereditary dependence. 

For ordinary differential equations, this method Is well understood 
by most people who are concerned either with the computational aspects or 
the qualitative theory of nonlinear oscillations. In the development of this 
method for retarded systems, the basic difficulty lies in the fact that 
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notions defined by the solutions of the equations cannot be described 
adequately in a finite dimensional space. The proper setting seems to be 
in an infinite dimensional space and, in the particular formulation given 
below, in a Banach space. To the author’s knowledge, Krasovski! [l] was 
the first to exploit such equations in this setting in the extension of 
Lyapunov's second method. 

The extension of the method of averaging to differential equations 
with retardation relies heavily upon the theory of linear equations with 
constant coefficients as developed by Shimanov [2, 3 ] and the author [4]. 

We will not give the details of the theory of linear systems, but merely 
apply the results to our problem, proceeding from its application to 
specific examples to the more general results. 

2. Notation. 

Let us digress for a moment and discuss the equation 
(2.1) i(t) ■ f(t, x(t), x(t - r)), r k 0, 

["•" represents the right hand derivative] from the point of view that is of 
interest to us in this paper. If x is in R n (the n- d imens i onal Euclidean 
space), f(t, x, y) is continuous in its arguments for all t, x, y, and 
<p is any continuous function mapping the interval [-r, 0] into R n , then 
for any t one can show that there is a function x(t Q , 9 ) which is defined 
on an interval [t Q - r, t Q + A) , A > 0, coincides with 9 on [t Q - r, t Q ] 
and satisfies (2.1) for t k t o « We call such a function a solution of (2.1) 
with initial value 9 at t Q . Furthermore, if f(t, x, y) is locally 
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Lipschitzian In x, y this solution is unique and depends continuously upon 
t Q and 9. If we find another solution with initial value i|r at t^ 
then the corresponding solutions might behave as in Fig. 1. 

The uniqueness property asserts that if two solutions coincide on 

any interval of length r, then they must coincide for all future time, but 

two distinct solutions may intersect many times on any interval of length r. 

This Hast remark su gg ests that the state of time t of a system described 

by (2.1) should be the collection of values of the solution on the interval 

[t - r, t], or the restriction of x to the interval [t - r, t], We 

designate this restriction by x^ (see Fig. 2). If we let C = C([-r, 0], R n ) 

be the space of continuous functions mapping the interval [-r, 0] into 

R n with the uniform topology, then a solution x(t Q , 9) of (2.1) yields 

for each fixed t £ t Q a mapping of C into C; namely, the mapping 

x t ( V *>• Trajectories of (2.1) are then defined as the collection of points 

(t, 9)) in R x C t Q S t < t Q + A, as indicated in Fig. 3. Hereafter, 

we assume solutions defined for all t g t - r; that is, A = + 00. 

o 9 

The above definition of trajectories of (2.1) yields a situation which 
is analogous to ordinary differential equations. However, the reader should 
realize that the situation here is more complicated. First of all, trajectories 
in general are only defined to the right of t Q and the mapping x^ is a 
smoothing operator if r > 0. In fact, for any t^ the mapping 
x^.(t o , 9) takes closed bounded subsets of C into compact subsets of C. This 
shows that x, (t , 9) cannot be a homeomorphism for r > 0 even if it is 
one-to-one* Secondly, the mapping x^(t^, qp) need not be one-to-one even vhen 


the uniqueness property holds. In fact, for the scalar equation 


*(t) = x(t - r) (l - x(t)), 

the solution x(0, <p) corresponding to an initial function <p with <p(0) = 1 
is such that x^(0, <p) = 1 for t fe r. Therefore, a subset of c([-r, 0], R) 
vhich is the translate of a subspace of codimension 1 is such that the 
corresponding trajectories all coincide after r units of time. 

If ve let F(t, <p) be a functional defined on [0, op) x C into R n , 
then a rather general hereditary functional-differential equation can be 
defined as 

(2.2) £(t) = F(t, x t ) 

where x^. is the restriction of x to the interval [t - r, t]. The 
discussion below is concerned with these more general equations, but we 
devote much of our time to more specific types. Equations (2.2) are 
certainly more general than (2.1) and certainly include (2.1) with the functional 
F defined by 


F(t, x t ) = f(t, x(t), x(t - r)). 

3. A convenient coordinate system. 

In this paper, we are interested in the oscillatory properties of 
perturbations of linear equations with constant coefficients. In ordinary 
differential equations experience has shown an understanding of oscillations 
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in perturbations of linear equations with constant coefficients is most 
easily accomplished by the introduction of a coordinate system which exhibits 
in an explicit manner the behavior of the unperturbed equation on the 
subspaces which correspond respectively to the eigenvalues which positive 
real parts, zero real parts and negative real parts. In this section, we 
indicate hew this same end can be accomplished for hereditary functional- 
differential equations. 

We will first discuss the procedure for the simple equation 

(3.1) u(t) = -au(t - r), <xr - ir/ 2 , r > 0, 

and the perturbed equation 

(3.2) £(t) = -ooc(t - r) + ef(t, x t ) 

vhere x^ denotes the restriction of x to the interval [t - r, t]. Later 
ve state a result for a more general equation* 

If a function e^ is a solution of (3*1), then X must satisfy 
the characteristic equation 

(3.3) X “ -c«~ Xr 

Since err = ir/ 2, it is not difficult to show that the roots of (3*3) all have 
negative real parts except for two which are equal to + ice. Furthermore^ 
every periodic solution of (3*1) must he of the form 


(3.4) 


a sin cet + h cos crt 


a 
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far same constant a and b and every solution of (3.1) is exponentially 
asymptotic to a function of the form (3.4) as t -» 00. 

Now let us interpret these remarks in the space C. Let u be a 
periodic solution of ( 3 * 1 ). Then for -r ^ 0 ^ 0 , 

\(e) = u(t + 0) = a sin a (t + 0) + b cos a (t + 0) 

= (a cos a t - b sin afc) sin a© + (a sin at + b cos at) cos a0 
d - f y^t)^ 0) + y 2 (t)<p 2 (0) 
where we have defined 

(3.5) ^(e) = sin O0, q>g(e) = cos aQ » -r S 0 S 0, 

and y , y^ are the corresponding coefficients of these functions. 

Now q>2 are linear independent elements of C and thus generate 

a two-dimensional linear subspace P of Cj that is, 

(3*6) P « (9 e C: 9 = aqp^ + txp 2 , a, b real). 

What the above computations have shown is that all of the periodic solutions of 
(3.1) must lie in P. Also, it can be shown that any solution of (3.1) approaches P 
exponentially as t -» 00. The paths in P are closed curves and the motion in 
time is described by y^t), y 2 (t). Notice that y^t), y g (t) satisfy the 
ordinary differential equations 



If we could find another subspace Q of C (which necessarily 


must be infinite dimensional) which is positively Invariant under the 
solutions of (3.1) and complementary to P in the sense that for every 
9 in C there exist unique elements q)p in P, <p^ in Q such that 
9 = q>p + 9^, then pictorially the motions in C would be as shewn in 
Pig. 4. 

The existence of such a space Q follows from the general theory of 
linear operators since P is an eigenspace of the semigroup of bounded linear 
operators u(t), t 0, defined on C by u(t)9 = u t (<p), where u(9) 
is the solution of (3*1) with initial function ip at 0. On the other hand, 
if Q can be described analytically, then we will be in a position to 
introduce a coordinate system in C which will provide a natural means of 
extending perturbation theory. 

This is accomplished by means of the equation 

(3.7) *<s) ■ ov( s + r ) 

"adjoint" to (3.1) with respect to the bilinear form 

0 

(3.8) (+, 9) = +(0)9(0) - a / +(| + r)9(|)d| 

-r 

defined for all + in C( [0, r], R) and 9 in C([-r, 0], R). This bilinear 
form has the property that if v is a solution of (3»7) defined for s 5 0 
and u is a solution of (3.1) defined for t S 0, then 

( v t* u t^ = cons ‘*' an ^ for ^ - 0. 



/ 
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If A is a column vector of dimension k vhose elements a^ 
belong to c([0, r], R 11 ) and B is a row vector of dimension m vhose 
elements belong to c( [-r, 0], R n ), then ve let (A, B) denote the 

kXm matrix vhose (i, element is given by (a^ b^). 

Equation ( 3 « 7 ) also has two linearly independent periodic solutions 
sin os, cos a s defined for s in (-«, co) . Define i| in 

C([0, r], R) by 


(3.9) 


and 


t x (e) - sin ad, t 2 (e) = cos aBf 0 ^ Q S t 


* c (t 1# i 2 ), * - 


u * = <P 2 ). 


Then a simple computation shows that the matrix 


(1, »> iSf ((t ± , 9j)) “ \ ( i ^ ) 


and it is nonsingular. It is convenient to define 


y" - Y(Y, 4>) _1 = Y 


H? 


7r 

2 


Since 


(**, •) = I, the identity. 


TT 

"2 
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We are new in a position to introduce a coordinate system in C and 
define the space Q complementary to P. In fact, for any 9 in C we let 

q> = *c +9, c * (**, <P)> 

which gives a unique decomposition of every el: 9 in C. The subspace Q 

complementary to P is defined by 

Q = {9 in C: (?*, 9) = 0). 

For the integral representation below, it is necessary to extend the definition of 
Q to piecewise continuous functions. It is clear* that this is possible and our 
decomposition is valid in this larger space. Hereafter, Q will denote this set. 

If x is a solution of (3.2) with initial value 9 at o, o' in (-.«,«), and 

(3.10) x t = *y(t) + x t , 9 * 4>b + 9, 

then y(t), must satisfy 

ftt) - By(t) + ef^O) f(4y(t) + x t ), y(a) - b , 

k.U) t 

) + € / u t _ T (x 0 ) f (®y( T ) + * x ) dT 


where 


- 10 - 



9, X o are in Q, u(t) is the solution of (3.1) with initial value t at 0, 
X D (e) » 0, -r s 0 s 0, X Q (0) = I, the identity. If 9 is in Q, then 
there are positive K, OC such that 

(5.12) KWH s Ke '°*ll»ll» 4 s °- 

This relationship expresses more precisely our stability property of the set 
P mentioned before* 

Now let us see what this coordinate is like for the general linear 
equation. 


i 

1 



0 

(3.15) u(t) = / [dT)(e)]u(t + 0), 

-r 

where ij is an n X n matrix whose elements are real functions of hounded 
variation, and the perturbed equation 

0 

£(t) 55 / [dT)(e)]x (t + 0) + ef(t, x t ). 

-r 


! 


(5.U) 
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We wish to Indicate how a coordinate system can be Introduced into C in such 
a way as to obtain a set of equations equivalent to (j.14) which is of the form 
( 3 . 11 ) with Uj.(q>) satisfying (3*12) and the corresponding matrix B having 
eigenvalues which coincide with the characteristic values of ( 3 « 13 ) which have 
real parts 2 0* The characteristic values of (3«13) are the roots of the 
equation 

(3.15) det [XI - /° Cdn(0)]e X0 ] = 0, 

-r 

and to any characteristic value, X, there is a solution of (3»13) the 
form e X ^b for some b and all t in (-«, «)* 

As Is to be suspected from the previous discussion, a basic role is 
played by the equation* 

0 T 

( 3 . 16 ) $(s) ■ - / [dTj (e)]v(s - 0) 

-r 

"adjoint" to ( 3 . 13 ) with respect to the bilinear form, 

0 0 

(3.17) (t,q>) c Y (0)<p(0) - / / t T (| - 0)[dit(0)]<p(|)d|, 

-r 0 

defined for all if in c( [0, r], R n ), <p in C([-r, 0], R n ). The 
characteristic values of the adjoint equation are the roots of the equation 

(3*18) det [XI - /°[dri T (0)]e X0 ] = 0, 

-r 


*If A is a matrix A T denotes the transpose of A. 
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and to each such root, X, there is a solution of (3.16) of the form 
e _XS h for seme constant vector b and all s in (-», 00). Notice that 
the solutions of equations (3.I5) and (3. 18) are the same. 


Suppose X , X^ are the characteristic values of (3*13) with 

real parts h 0. There are only a finite number, say m, of linearly 

independent solutions of (3.I3) of the form Pj(t)e^J^ where the 

are polynomials. Let 4 = (9^...,^) where 9 .^ , 9^ are the 

restrictions of these functions to [-r, 0], Similarly, there are only m 

_k 

linearly independent solutions of (3.16) of the form lj(s)e J where 

the q^ are polynomials. Let ¥ = col where +1# • • • 1 l f m are 

the restriction of these functions to [0, r]. 


It follows directly from the differential equations that there is a 
square matrix B vith only the eigenvalues X 1 ,«..,X fc such that 

(3.19) *(0) = *( 0 )e m , -r * 0 ^ 0. 

Furthermore, one can show (see Hale [4]) that the matrix (y, 4) is 
pnpnipgiilaT and, therefore, by a change of the basis Y, one can take 
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(*, *) to be the identity. Finally, the transformation (3.10) with 

f, ¥ as above applied to (3»l4) yields an equation of the form ( 3 «ll) with 

^( 9 ) satisfying ( 3 » 12) and the matrix B given by (3.19)» 

V. Perturbation theory. 

It was indicated in the previous section that there is a transformation 
x t = 4 >y(t) + x t , x t in Q, which takes the general system (3.14) into an 
equivalent system of the form 


(*.l) 


${t) “ By(t) + e* (0)f(t, 4y(t) + x t ), y(<r) = b 

t 

\ m + 6 / U t . T (X 0 )f(T, ®y(T) + X T )dT 

o 


where 9 * #b + 9 , 9 in Q, the eigenvalues of B have nonnegative real 
parts, u(9) is the solution of (3*13) with initial values 9 at 0 and 


|lu t (9)|| £ Ke" afc ||9||, t £ 0, K > 0, a > 0, 


for any 9 in Q. 

Equations (4.1) are now in a form which is very similar to that which 
is encountered in the theory of oscillations in ordinary differential equations. 
One can show that any solution of (4.1) which is bounded on (- 00 , ») must be 
of such a nature that x^. = 0( c) as e -> 0. Consequently, if our analysis is 
based upon an approximation procedure which can be Justified to be correct by 
investigating only the terms of order €, then the basic problem lies in the 
investigation of the ordinary differential equation 
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(4.2) HO = Sy(t) + eir^oJfCt, *y(t)). 

The analysis of (4.2) is well understood and usually proceeds by the in- 
troduction of convenient combinations of polar coordinates and rectangular 
coordinates and the application of averaging procedures and successive 
approocimat ions . 

For simplicity, let us mate the assumption that all the characteris- 
tic values of (3*13) have nonpositive real parts. Then B in (4.1) has 
all eigenvalues purely imaginary and a combination of polar and rectangular 
changes of coordinates in the components of y (see the examples in 
section 3 for the types of coordinates involved) leads to a set of equations 
of form 


£ ■ d + e0(t, 5, p, x t ) 

(4.3) f ~ € R(t, £, Pi x -^) 

t 

*t = u t-o- ( ^ +€ ^ u t-T (f o )P(T ' k p ' *r )dT 

O' 

where £ is a p -dimens ion vector, p is a q-dimensional vector, x^. 
is an element of the Banach space C. The vector d is a constant vector 
with positive components and the functions 8 , R, F are multiply periodic 
in the vector £. 

Assume that the functions 8 , R, F with arguments t, 5 , p, 9 have 
continuous second derivatives with respect to £, p, <p and are almost 
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periodic in t uniformly with respect to £, p, 9 in some set. Let 


= + ..., Cp + T ) 


and assume that 


t T def 

(kA) lim = / R(t + t, 5 + p,o)dx = R (p) 

T -» <» 0 0 


is independent of t, £. We define the averaged equations associated with 
(4.3) to he the equations 

(4.3) £ « € R c (p). 

Hotice that the averaged equations (4.5) are obtained from R(t, £, p, 0) 
and, therefore in a specific problem, they arise from an investigation of 
the ordinary differential equation ( 4 . 2 ). 


Theorem 4 .1. If system (4.3) satisfies the conditions enumerated above and 
if there exists a vector p Q such that ® 0 (p 0 ) = 0 and the eigenvalues of 
the matrix ^® 0 (P 0 )/^P have nonzero real parts, then there exists an e Q > 0 
and functions g(t, 5 , e), h(t, £, e), 0 ^ g in R q , h in C, 

g(t, 5 , 0 ) ® p o , h(t, £, 0 ) = 0 , multiply periodic in 5 and almost 
periodic in t such that the set S £ , 0 S e £ e^, defined by 

s € ((t, p, 9): p - e(-fc, 5 ,< 0 , 9- *(t,5,e), 

-» < t < », - 4n < < "> d® 1 * 2 * •••* p) 
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is an integral manifold of system (4.3). If the functions 0, R, F 
are independent of t (or periodic in t of period ®), then the func- 
tions g, h are Independent of t or periodic in t of period <o). 
Furthermore, if all eigenvalues of dR Q (p o )/c)p have negative real parts, 

S £ is asymptotically stable for 0 < e £ e Q and if one eigenvalue has 
a positive real part then S £ is unstable for 0 < e £ e Q . 

We merely give an indication of the proof of this theorem since it 
is so analogous to the proof for the case or ordinary differential equa- 
tions given in Bogoluibov and Mitropolski [ 5 ] and Bale [ 6 ]. 

In [ 6 , Ch. 12], it is shown that there is a function v(t, £, p, e), 
multiply periodic in £ and almost periodic in t such that the trans- 
formation 

P -»P + ev(t, t , P> € ) 

applied to the equation £> = e R(t, £, p, 0 ) yields a new equation of the 
form !> « e R c (p) + e Ri(t, £, p, e ) where R^(t, £, p, € ) is zero for 
e = 0. Consequently, if this transformation is applied to (4. 3)> we obtain 
a system of the form 

i - d + e ©j(t, £, p, x t , e ) 

f> - € R c (p) + € R x (t,5,p,e) + e R 2 (t,£,p,x t ,e) 

t 

X t = + € / li t ^(X o )F 1 (T,5,p,X T ,€)dT 

u 
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vbere 8^, F 1 are the sane types of functions as 6, F and 

£> P# 0) ■ 0, Rg(t, 5>P>0^e) = 0, |Rg(t, <^e ) - Rg(t, £,p,i|r,€ ) { - 
K||9 - i(r|| for some constant K and p, f Isa hounded set. 

One nov proceeds in a manner completely analogous to that given in 
[5\, [6] to show that the functions g, h mentioned in the theorem are the 
fixed points of an integral operator. The stability of the integral mani- 
fold must he investigated separately and is easily supplied using the 
ideas developed in [8] in connection with a saddle point for functional- 
differential equations. 

As in [ 5 ], [6], one can also prove 

Theorem 4.2. Suppose the averaged equations (4.5) have a nonconstant 

periodic solution p = p°(t) of period T such that q-1 of the characteris 

tic exponents of the associated linear variational equations have nonzero 

real parts. Then there exists an e Q > 0 and functions g(t, £, +, e), 

h(t, t,, if, «), 0 £ e - e Q > g in R q , h in C, g(t, £, +, 0) = 

p°(+), 0 ^ f £ T, h(t, £, if, 0) = 0, multiply periodic in £, periodic 

in + of period T and almost periodic in t such that the set S £ , 

0 £ e £ e defined hy 
o 

S e *= ((t,£,p, 9 ): p = g(t ,£,+,€), 9 = h(t,£,+,e), 

.»< t <»;-«< i ■ h 2 > •••> PJ 0 s t s 1 ) 

is an integral manifold of system (4.3). If ©, R, F are independent of 
t (or periodic in t of period 00 ), then the functions g, h are 
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inde pendent of t (or periodic in t of period oj). Furthermore, the 
stability properties of S £ are the same as those of the periodic 
solution p°(t) of (4.5) • 

We now state some important corollaries of these theorems before 
turning to specific examples. Consider the equation 

' .6) £(t) = €f(t, Xj.) 

where e > 0 is a parameter, f (t, 9) is almost periodic in t uniformly 
with respect to 9 in some subset of C([-r, 0], R n ), and has a con- 
tinuous second Frechet derivative with respect to 9. Let 

(4.7) f G (q>) = U* £/*(*, 9)^. 

Hal arm y [7] has discussed, for small e, some of the relationships between 
the solutions [on the interval (0, <*03, of x(t) = €f Q (x^.) and the solu- 
tions of (4.6) for the case in which the retardation interval is of order 
€. We now show that Theorems 4.1 and 4.2 imply that his results and even 
more are valid without any restriction on the retardation interval. In 
fact, we can prove the following two theorems. In the statement of these 
theorems, y sometimes denotes a vector in n -dimensional Euclidean space 
and sometimes a vector of constant functions in C([-r, 0], R n ) , but it 
is clear from the context which meaning is implied. The averaged equations 
of (4.6) are then defined to be the ordinary differential equation 
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(4.8) f «= €f o (y). 

Theorem 4.3. If the averaged equations (4.8) have an equilibrium 
point y Q such that the matrix of coefficients of the linear variational 
equations has no eigenvalues on the imaginary axis, then, for e suffi- 
ciently small, (4.6) has a unique almost periodic solution x = g(t, e) 
in a neighborhood of x = y Q , g(t, 0) = y Q , and the stability properties 
of g are the same as the stability properties of y Q . 

Theorem 4.4. If (4.8) has a nonconstant periodic solution 
y — y(°^(t) of period T, such that the linear variation equation has 
n-1 of its characteristic exponents not on the imaginary axis, then, for 
e sufficiently small, there exists a function g(t, £, e) in C, almost 
periodic in t uniformly with respect to £, periodic in £ of period 
T, 6 0** £> 0) = y(°\ y^ Q \e)= y(°' J (£ + 9), -r 5 0 5 0, such that the 

surface S £ in C X (-*», 00 ) defined by 

S £ = {(<p, t): qp = g(t, 5, €), 0 s 5 S I, < t < •») 

is an integral manifold of (4.6). Furthermore, S £ is unique in a neigh- 
borhood of S Q = {(<p, t): qp = y^ 0 ^, 0 ^ g T, .<* < t < «) and has 

the same stability properties as S Q . If f in (4.6) is independent of 
t, then g is independent of t, and if f is periodic in t, then 
g is periodic in t with the same period. 
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To show that these results axe consequences of Theorems 4.1 and 
4.2, we proceed as follows. For any 9 in C([-r, 0], R n ), the de- 
composition 

9 = h + 9> 

is unique if h is the constant function vhose value is 9 ( 0 ). If, in 

(*. 6 ), 

= y(t) + x t 

then 

t(t) - ef(t, y(t) + x t ) 
t 

+ € 1 \-T ( *o )f(T ' y(T) + *r )dT 

0 

and 

I^OSII SKa-“ t ||?||, «>0 

for some positive K, a. This last relation is obviously true in this case 
since 1 ^( 9 ) = 0 for all 9 and t S r. SyBtem (4.9) is a, special case 
of (4.?) and one' obtains Theorems 4. 3 and 4.4 from Theorems 4.1 and 4.2. 

5 . Some specific examples . 

Let us first discuss the oscillatory properties of equation (3*2), 


namely, the equation 
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i(t) *= -ooc(t - r) + ef(t, x t ) 
or = v/2, € > 0. 


life have seen in section 3 that this equation is equivalent to the 
system (3.1l), that is, the equation 

= -ay 2 + eir|i 2 f(t, *y + x t ) 

*2 - °yi + *y + ^fc) 

(5.2) ^ = H^?) + € + *r> dT 

*(0) = (<P,.(0), <fc(e)) ^(sin as, cos a0), -r s e S 0 



If ve let 


(5.3) 


y 1 = p sin at, 
y 2 *= -p cos a£ 


then system (5.2) hecomes 

j«:l + ^-(irf cosat; + 2f sin a£) d i f 1 + i, P> x t ) 

(5.^) £ ■ ep 2 Grf B i n a£ - 2f cos a£) e R (t > £* 

x x . ^ci) + £ /V_A )f(T ' «y( T > + *,> 4T 
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vhere f = f(t, *y + x^.) and y is given in ( 5 . 3 ). 

If the function f is almost periodic in t uniformly with res- 
pect to the other arguments, then the averaged equations for this parti- 
cular case are 


A *= € R Q (p) 

R (p) = lim i / \ifa f sin a(£ + t) - 2f cos a(£ + T)]dr 
T -*«* 0 

(5.5) f = f(t + t, ptq^sin a(£ + t) - <p 2 cos a(£ + *)]) 

9 ^( 0 ) = sin ad, <^( 0 ) = cos a0, - r £ 0 £ 0, 

where we always assume that this limit is independent of t, £. 

notice that equation (5«5) are the same equations that are obtained 
by introducing the polar coordinate transformation (5.3) to the ordinary 
differential equation 

■ - ay 2 + ^ 

(5.6) 2 , 

* Wx + 2eu *(*> *ar) 

and then taking the average. 

Let us now tr.:- - sane specific functions f in (5-1) to show that 
important informs, . rs obtained by this method. 

Example 5.l(Pinney 19 J f - - ipc(t-l) + Bx^(t-l), q = tt/ 2, r = 1. 

Since x(t-l) = *(-l)y(t) = -y^t) = - sin(ir/2), it is easy to 
check that the average R Q (p) in (5*5) is 
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R 0 (p) = m>(l - ^ ) , 

■where K is a positive constant, and the averaged equations (5.5) are 

£ * € *5fP(l 

If x& > 0 the averaged equation has an equilibrium point p Q = */4y/3P 
which is asymptotically stable if x > 0 and unstable if y < 0. Con- 
sequently, for e > 0 and sufficiently small. Theorem 4.1 asserts for 

yp > 0 the existence of functions g (£, e), h(£, e) (these functions 

* . 

are . impendent of t since f is)^ g(£, 0) = *1 4y/3P, h(£, 0) « 0, 

per in £ of period 2rr/a such that S £ = {(t, £, p, 9): 

P - e), 9 = h( £, e), -o»<t<«, OS £ S 2n/a} is an integral 

manifold of (5.4) which is stable for y < 0. From (3.10) and (5* 3) , this 
implies that"^s((t, 9) : 9 = (j^sin a£ - <j^ cos ec£)g(£, €) + h(£, e), 

0 £ £ S 2jr/a) is an integral manifold of our original system. Such a 
cylinder T £ in R x C obviously corresponds to a periodic solution of our 
Bystem which is stable if y > 0 and unstable if y < 0 and has an ampli- 
tude approximately equal to ^4y/3P« The approximate period a> is obtained 
by solving the equation 

£ » 1 + eu 2 (r sin tt£/ 2 - fr 0 5 sin 3 ir£/2), P Q = */ 4y/3p, 
and determining . o> so that £(t + oo) = $(t) + 2ir/a. 
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It is interesting to note that the second order system (5*6) for 
this example is actually equivalent to a second order scalar differential 
equation. The method of averaging should then allow one to obtain an 
"equivalent" linear second order equation in the sense of Krylov-Bogoliubov . 
This should in turn lead to methods which will yield important information 
about equations with retardation when € is not snail — describing func- 
tions, etc. So far, this has not been exploited. 

Example 5.2. Consider the equation 

(5.7) £(t) c - [| + €q(t)]x(t - 1)(1 - ex 2 (t)) 

where tj is almost periodic in t. 

This is a special case of ( 5 . 1 ) with a = 7T/£, r = 1 , and 

f = | x 2 (t)x(t - 1) - q(t)x(t - 1) + eri(t)x 2 (t)x(t - 1). 

Using the fact that 

x(t) = *(0)y(t) ■ y g (t) = - p cos (irt/2) 
x(t-l) » *(-l)y(t) = - y x (t) = - p sin(uV2) 


the averaged equation (5*5) becomes 
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T 

fi * - g p 2 ), = lim | / rj(t)dt 


provided 


kj n(t) 
r -too 1 o 


( cos *-(£ + t) 
cos Tr(t + t) 


dt = 0. 


Consequently, if > 0 and € > 0 is sufficiently small there exists 
a stable integral m a nif old of solutions of (5.7) whose parametric repre- 
sentation in C is a lm ost periodic in t and periodic in £ and for 
e « 0 is given by 


p sin vt/2 


L -p ft cos TTt/2 



If tj(t) is independent of t, then the parametric representation 
of the integral manifold is independent of t and one obtains a nonconstant 
periodic solution of (5. 7) with amplitude approximately '/8q/V. Jones [IX)] 
has discussed periodic solutions of (5*5) with tj independent of t and 
even more general equations. 


Example 5 . 5 * Consider the system 


(5.8) 


*l(t) - x 2 (t) 

i 2 (t) = + e[l - Xj^ 2 (t - r)]x 2 (t) 
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where a > 0 , r £ 0 , e > 0 are parameters. For r = 0, this is 
van der Pol*s equation. By using the preceding theory, we will invest! 
gate the existence and stability of limit cycles of (5*8) for e small 
If <f^ are the vectors in C([-r, 0], R ) defined hy 


cos ad 


fr sin a$ 


9,(0) = / \ , q^(e) 4° J, -rseso 

,-a sin a QJ Vcos a© 


and 4> = (<^> 9^) then the transformation 


x t = *y(t) + x t , y(t) = x(t), y =j ], x = 


applied to (5.8) yields the equivalent system 


(5.9) + e [1 - (y x COB ar - i y 2 sin ar + x t (-r)) 2 ]y 2 (t) 

t _ 

x t * ^(9) + e / \_ T (X Q )F(y(T), ^)dT 


where 1^(9) has the same meaning as in the previous sections and F 
two-vector whose specific form is of no particular interest here. 

If we introduce the polar coordinates 


y x « p sin a£ 
y 2 «= pa cos a£ 


is a 
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Into (5*9) and set x^ = 0, we obtain 

C » 1 - ^ sin 2a£[l - p 2 {sin a(£ - r)} 2 ] 

0 = €p cos 2 a£[l - p 2 (sin a(£ - r)} 2 ] 

The averaged equation is 

& - [1 - |(1 - | cos 2ar)p 2 ]. 

This equation has an equilibrium, point p Q = y = 1 - (cos 2ccr)/2 > 0 

for every value of r and the linear variational equation relative to p Q 
is £ = - 2ep. Consequently, Theorem 4.1 implies as in Example 5*1 the 
existence of a stable periodic solution of (5.8) with amplitude approxi- 
mately */2/V. 

For physical examples of retarded equations and the importance of 
oscillatory phenomena, see Chapter 21 of N. Minorsky [11]. 
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